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^ ■ ABSTRACT 

The presence of an imposed external magnetic field may drastically influence the 
■ structure of thin accretion discs. The magnetic field energy is here assumed to be in 

psj ' balance with the thermal energy of the accretion flow. The vertical magnetic field, its 

toroidal component 5'°'' at the disc surface (due to different rotation rates between 
disc and its magnetosphere) , the turbulent magnetic Prandtl number and the viscosity- 
' alpha are the key parameters of our model. Inside the corotation radius for rather small 

5*°'' the resulting inclination angle i of the magnetic field lines to the disc surface 
normal can exceed the critical value 30° (required to launch cold jets) even for small 
magnetic Prandtl numbers of order unity. The self-consistent consideration of both 
magnetic field and accretion flow demonstrates a weak dependence of the inclination 
("dragging") angle on the magnetic Prandtl number for given surface density but a 
strong dependence on the toroidal field component at the disc surface. 

A magnetic disc is thicker than a nonmagnetic one for typical parameter values. 
The accretion rate can be strongly amplified by large B^°^ and small magnetic Prandtl 
number. On the other hand, for given accretion rate the magnetised disc is less massive 
Q ■ than the standard-alpha disc. The surface values of the toroidal magnetic fields which 

| H \ are necessary to induce considerably high values for the inclination angle are much 

C/3 . smaller than expected and are of order 10^'^ of the imposed vertical field. As the 

' innermost part of the disc produces the largest the largest radial inclination can 

be expected also there. The idea is therefore supported that the cold jets are launched 
. . only in the central disc area. 

rS ' 
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1 INTRODUCTION 

Accretion flows is assumed to occur in a variety of astro- 
physical objects. The structure of magnetic accretion discs 
is the problem of major importance, for example, for under- 
standing of the origin of astrophysical jets. If the inclina- 
tion, i, of the field lines from the vertical exceeds 30°, the 
plasma can be accelerated when spiralling along the field 
line (Blandford & Payne 1982; Lynden-Bell 1996; Camp- 
bell 1997; Krasnopolsky, Li & Blandford 1999). Accretion 
disc-driven magnetocentrifugal winds have been widely used 
to model the astrophysical jets (e.g. Ustyugova et al. 1999, 
Ouyed & Pudritz 1999, Krasnopolsky, Li & Blandford 1999). 

The present study is motivated by a series of papers 
dealing with the interaction of accretion discs with external 
magnetic fields (Livio & Pringle 1992; Lubow, Papaloizou 
& Pringle 1994; Bardou & Heyvaerts 1996; Reyes-Ruiz & 



Stepinski 1996; Ogilvie 1997; Ogilvie & Livio 1998; Campbell 
1998; Campbell & Heptinstall 1998a, 1998b; Brandenburg 
& Campbell 1998). In the presence of a vertical magnetic 
field, the vertical gradient of angular velocity generates the 
toroidal field due to the stretching effect. This field may 
generally infiuence the vertical and radial structure of a disc 
and alter the angular momentum transport. If the magnetic 
field is sufficiently strong it even will infiuence the rotation 
law (Ogilvie 1997; Ogilvie & Livio 1998). 

Our previous paper (Shalybkov & Rudiger 2000, here- 
after paper I)) considered the self-consistent steady-state 
structure of the polytropic accretion disc in the presence of 
a vertical magnetic field. In the present paper we will refuse 
from simplified polytropic equation of state and solve the 
full system of MHD equations for the vertical accretion disc 
structure. So better comparisons with the theory of the ver- 
tical structure of accretion discs without magnetic field are 
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possible. As known, such computations, for given viscosity- 
alpha and given opacity law, yield the accretion rate M for 
any possible column density E. The same is done here for 
an accretion disc threaded by a magnetic field with a given 
vertical component B^^'^^ and a given toroidal component 
-8*°"^ which, of course, have opposite signs if we are within 
the corotation radius*. Such a configuration can only exist 
for a certain accretion rate and with a certain distribution of 
a radial magnetic field component Br. The value of the lat- 
ter taken at the surface defines the inclination angle i known 
from the jet theory. Along this way the well-known dragging 
problem - to find the radial inclination i - has been unified 
with the theory of the vertical structure of accretion discs. 

There is an estimate for the toroidal component _B*°' 
due to the different rotation of halo and disc. If a disc halo 
with a high conductivity exists then the surface value of B^ 
will become large. After Campbell (1992) the shear between 
the rigidly rotating halo and the accretion disc induces a 
toroidal magnetic field of 
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with f2, as the stellar rotation rate. The numerical integra- 
tion of the induction equation for a uniform and a dipolar 
field and a plasma halo (with conductivity of 10 times of 
the disc conductivity) confirms this result with 7 ~ 1 (Elst- 
ner & Riidiger 2000). The toroidal surface field changes its 
sign at the corotation radius where f^Kep equals the stellar 
rotation rate fi*. The magnetic torque results as negative 
inside the corotation radius and positive outside the corota- 
tion radius. For a disc embedded in vacuum, of course, the 
7 in (Q) vanishes. We consider the 7 as representing the un- 
known halo conductivity. We shall show here that we only 
need very small 7 in order to produce inclination angels i of 
the interesting value of 30° and more - independent of the 
the magnetic Prandtl number Pm. We also can take from 
(^) and from the simulations by Elstner & Riidiger that the 
ratio 
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- which is negative (positive) inside (outside) the corota- 
tion radius grows (by absolute value) inwards. This fact will 
play an important role in the philosophy of the presented 
research. 



2 BASIC EQUATIONS 

Consider the structure of an axisymmetric disc in a steady- 
state regime; R, (j>, and z are the cylindrical coordinates. 
The resulting kinetic and magnetic equations are given in 
paper I as the Eqs. (1)...(13) there. Here we have only to add 
the energy equation. All generated energy is assumed to be 
transfered from the disc by radiation neglecting convection, 
and the disc is optically thick. Then the energy equation is 
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* At the corotation radius the component 5'°'' vanishes; for a 
nonrotating central object the corotation radius is in the infinity 
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where the energy flux components are 

IGctT^ dT _ IGctT^ dT 
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Skp 9_R' 3«;p dz 

with (7 is the Stefan-Boltzmann constant, T is the tempera- 
ture and K the Roseland mean opacity. 

The equations must be supplemented by relations spec- 
ifying the gravitational potential, the equation of state, the 
opacity, the viscosity, and the magnetic diffusivity. We ne- 
glect the self-gravitation of the disc so that 

(i?2+ ^2)1/2' 

where G is the gravity constant and M, is the mass of the 
central object. We adopt here the ideal-gas equation of state 
P = TZpT / fi (TZ is the molar gas constant and p is the mean 
molecular mass) and the opacity may fulfil a power law, 
K — kop'^T^ , with the constant quantities kg, 7 and 5. This 
includes the cases of Thomson scattering opacity (fco ~ 0.4, 
7 = 5 = 0) and Kramers opacity (fco ~ 6.6 • 10^^, 7 = 1, 
5 — —3.5). A Shakura-Sunyaev parameterisation is used for 
the turbulent viscosity, 



pi^T = ass 



(6) 



where ass is a constant (Shakura & Sunyaev 1973). For the 
magnetic diffusivity we assume that the magnetic Prandtl 
number Pm = vx/rjT is constant. 

FoUowing Regev (1983) and Kluzniak & Kita (2000) 
we scale all quantities by their correspondent characteristic 
values. This will make the equations dimensionless and al- 
lows to compare the relative significances of each term. The 
radial distances are scaled by some characteristic radius, R, 
and vertical distances by a typical vertical height of the disc, 
H. We represent the angular velocity in units of the Kep- 
lerian velocity at the characteristic radius, = GMt./R?' . 
Using these three characteristic quantities all others can be 
defined. The typical sound speed is Cs = J/fi, the typical 
viscosity and magnetic diffusivity v = fj — CsH , the typical 
temperature is T = fi/TZc^, the typical pressure is P = pcf, 
the typical energy flux is _F = pel, and the typical magnetic 
fleld is _B = {fiopy^'^Cs. The last three quantities use typi- 
cal density, p, which transforms the equation for the vertical 
energy flux to the dimensionless form, i.e. 
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Note that the typical magnetic field deflned above B gives 
Va = B/{popY^^ = Cs, where Va is the typical Alfven ve- 
locity. The last relation defines the magnetic field energy as 
in balance with the accretion flow energy. The dimension- 
less equation of state takes the simple form P = pT, where 
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the same symbols as before are used for the dimensionless 
quantities. 

Let us define the parameter e 



e = H/R = Va/^R. 



(8) 



As in paper I we consider a geometrically thin disc (e <C 1) 
and expand all variables in power of e, i.e. 

A{R,z) ^ Ao{R,z) + eAx{R,z) + .... 

Using the same symbols for the normalized quantities, the 
full system at the leading order of e takes the form 
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where B^o = 5^°"^' and f^o = R~^^^ is the normalized Kep- 
lerian velocity. Equations (L)....(n3^ are identic with those of 
paper I. Note that after (hi) and(|l5[) all energy is assumed 
to be transported in the vertical direction. 

The system (|)....(^ depends on R only as parameter 
and we can use the real radius R as scaling constant R, the 
real disc half-thickness H as scaling con stant H a nd the real 
Keplerian angular velocity f^Kcp = \/ GKh / R^ as scaling 
constant Q,. The system takes thus the final form 
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and 



F ■ 



(22) 



if also the suffices and 1 are dropped. These relations form 
a nonlinear set of differential equations for P, T, F, ur, u^, 
Br, B^ with ass, Pm, and B^°'^^ as parameters. There are 
also two implicit parameters, i.e. R and the central mass 
M*. Solving the equations for given 7? and M, the verti- 
cal structure for given ass, Pm and B^"^ is found with 10 
reasonable boundary conditions. 

For discs symmetric with respect to the midplane {z — 
0) the physical quantities such as fl, ur, T are even functions 
of z while Br, B^, F are odd functions of z. Such symmetries 
provide the boundary conditions at the midplane {z — 0) as 



Br — B^ — duR/dz = du^/dz — F = 0. 



(23) 



At the disc surface, 2 = 1, it no external pressure is al- 
lowed. The simplest boundary condition for the tempera- 
ture, T(l) — 0, is used. It refiects the fact that for optically 
thick disc the surface temperature must be much smaller 
than the temperature at the disc midplane. Then the first 
of the above equations immediately gives ur{1) — for 
T(l) = 0. We additionally fix the toroidal and vertical mag- 
netic field at the surface as the last two boundary conditions. 
It means that the angle /3 is fixed and also the accretion rate 
(see (Eih and (|35|)). Hence, the boundary conditions at 2 = 1 
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where 5*°'' and i3^<=''* are free constants. As usual in the 
theory of the vertical structure of accretion discs, the ac- 
cretion rate (or the column density E) remain the only free 
parameter. 



3 PHYSICAL PARAMETERS 

It is useful to connect our normalized quantities with the 
real physical values. The real physical parameters of the 
problem are the distance to the central object, R, its mass 
Mt, the accretion rate M and the inclination angle i of the 
magnetic field lines to the rotation axis. Note that we can 
choose any of the accretion rate, the surface density and the 
disc half-thickness as the free parameter of the problem. All 
these quantities are connected with each other (see below). 
Accretion rate is more convenient as a parameter from the 
physical point of view, but the disc half-thickness is more 
natural for numerical calculations. 

The accretion rate Af and the surface density E as func- 
tions of half-thickness H are 

M = AirRpc^HKl, (25) 

where the p is given by (|^) and the normalized accretion 
rate M is 



M = — URopodz 
Jo 

and 



E = 2p//E, 



(26) 



(27) 
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where the normahzed surface density is 

E = /" podz. (28) 
Jo 

The disc thickness and the surface density as functions of 
M are 

i-5 / ,v \ 7+2\ 10-24 + 37 

where the e after (H) is 



and the surface density 
E = 



7+2\ 8-24 + 2-y 

(30) 
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Half- thickness and accretion rate as functions of E are 

1 

E \ 2+7 
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Integrating Eq. (119) over z one finds 
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so that ( |2i5| ) turns into 
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The accretion is thus provided by the surface toroidal mag- 
netic field in our model in strong contrast to the standard- 
alpha disc theory. According to the latter, the surface den- 
sity Eo of the disc is Eo = M jZixv with v = 2assCBHo and 
Cs = 2HonKep, where M and -ffo are the accretion rate and 
the disc half-thickness without a magnetic field. 1^ The Eo 
and Hq as functions Af*, R and M are 



8-24\ 10-24+^7 
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We are here interested to know the magnetic influence on 
the accretion rate, the disc height and the surface density. 
We shall consider two ratios for each quantities, for given 
accretion rate or for given surface density, resp.. Neverthe- 
less, all these characteristics are related to each other, and 
consideration of one ratio is sufficient. For given accretion 
rate, it results 
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and for given surface density 
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Our calculations will reveal the latter quantity always ex- 
ceeding unity so that for given column density E the related 
accretion rate M is higher for magnetised disc; and for given 
accretion rate the related column density is smaller than 
without magnetic field. 

The Mach number of the radial fiow is 



Ma: 



Cs 



yl/2 ' 

where we have used the ideal gas sound speed = 5/3 P/p- 
The inclination angle i of the radial magnetic field com- 
ponent at the surface Brs to the rotation axis is 



tani : 



(43) 



Its determination for prescribed accretion rate is the central 
point of the present paper, the often formulated dragging 
problem corresponds to the theory of the vertical structure 
of accretion discs in external magnetic fields. 

Below we shall use the constants R — lO^^cm, M* = 
IMq, p = 0.6, fco = 6.6 ■ 10^^ 7 = 1, 5 = -7/2. The Kepler 
velocity is f^Kep = 1.2-10~^s~^. The normalisation constants 
with e = 0.01 as a typical value^ are p = 2.8 ■ 10~* g/cm^, 
Cs = 1.1-10^ cm/s, f = 8.7-10^ K,i> = fj = 1.1 lO" cmVs, 
F = 3.7 ■ 10^'^ erg cm~^s^^ and the unit for the magnetic 
field is B = 590 G. 
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4 RESULTS AND DISCUSSION 

The systems of ordinary differential equations (|l^). . .(|22[) 
with boundary conditions (^3|) and (|2^) are solved by a re- 
laxation method (see Press et al. 1992). The dependent vari- 
ables are P, T, F, ur, u^, Br, and B^. The parameters are 
ass, Pm, M, iJ™"^' and B*""^. Only Kramers opacity is used. 



t We shall always neglect the factor 1 — (R^/R)^^^ with /?« 
the central object radius. 



t It can be easy calculated according to ^ or (g) and models 
values specified above and in Table 1. 
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Figure 1. The vertical disc structure for ass = 0.1, Pm = 1 and 
^vert ^ I rpj^g /3 is -4 ■ 10-4 (dot-dashed), -2.6 ■ 10-=* (short- 
dashod), -8.2 ■ lO'^ (dotted), -1.9 ■ 10"^ (solid), -3.5 ■ lO'^ 
(dot-dot-dot-dashed), —5.1 ■ 10"^ (long-dashed). The values of 
M, S and i for each curve are given in Table hi 



The calculation demonstrated that for a fixed -6*°"^ one 
always can adjust -B-Rs to fulfil the condition (^). On the 
other hand, for fixed Brs there is some critical value of the 
ass for given Pm, B^"^ and fixed by (H). It was pos- 
sible to find the solution only for qss < ctcr- The Ocr ~ 0.4 
for Pm = 1, -B""* = 1 and Brs = 1. Due to the nonlin- 
ear character of the system the reasons for this numerical 
problems are difficult to understand. Nevertheless, the real 
physical interest requires calculations with ass < 1- We have 
restricted ourselves to calculations for ass < 0.1 only. 

The vertical disc structure is illustrated in Fies. |l|and|[ 
We performed the calculations for vertical disc structure as 
function of Pm and for one of three parameters {M, E 

and i) fixed. Two values of viscosity-alpha, i.e. ass = 0.01 
and ass=0.1, are used always. The qualitative behaviour 
of the solutions is very similar for both cases. The profiles 
m Figs. and | are labelled by the toroidal inclination /3. 
We could them also label either with the nondimensional 



Figure 2. The same as on FigJll, but for agg=0.01 



accretion rate M or the nondimensional surface density E 
according to Table 1. 

Almost the entire variation of the disc variables hap- 
pens for < z < 0.5 hence always we have 'atmospheres' for 
z > 0.5. The atmosphere becomes thinner for increasing val- 
ues of the vertical field. The flow is subsonic (Ma < 1) and 
the applied surface toroidal magnetic field is much smaller 
in comparison with iJ""' (|^| < 0.1) for all calculated cases. 
The density and the temperature weakly depend on the tur- 
bulence parameters ass and Pm for large ass- The accretion 
fiow at the disc midplane can be comparable with the sound 
velocity. However, due to the weak dependence of the tem- 
perature on ass, the Mach number behaves almost linear 
with ass and therefore Ma ^ 1 for ass^ 1- 



4.1 Accretion rate and surface density 

The accretion rate for magnetic discs drastically increases 
in comparison with nonmagnetic discs (see Figs. ^ and 
Nevertheless, the accretion rate for magnetic discs can be 
even smaller than for nonmagnetic discs for large Pm and 
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Table 1. The model parameters (Pm= 1, B"" 
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Figure 3. The magnetic-disc accretion rates (normalized with 
accretion rates without magnetic fields) as a function of a) the 
vertical field S™"^' for ass=0.01, /3 = -2 - lO'^ (solid) and S = 
1.4 g/cm^ (dotted); b) the inclination angle i for B^'^'^ = 1, 
ass=0.1 (solid), ass=0.01 (dotted). ^ is defined by (^. Pm=l 



small i (the last case is the same as the case with large B^"^ 
when E fixed). 

Due to the direct relation between accretion rate and 
surface density (Fig. ^ the dimensionless surface density 
dependence on parameters is the same as for the dimension- 
less accretion rate. The surface density for magnetic discs 
is usually smaller than for nonmagnetic discs for the same 
accretion rate. The magnetic disc radial velocity is greater 
than nonmagnetic disc radial velocity (the accretion rate is 
greater in magnetic disc). The surface density does not de- 
pend on ass because the toroidal magnetic field is always 
small and we can neglect its influence on the pressure. 



4.2 Disc thickness 

In Fig. 1^ the thickness for magnetic discs is shown in compar- 
ison to the nonmagnetic case. The magnetic disc is thicker 
than the nonmagnetic one for all calculated models. Never- 



Figure 4. The accretion rate as a function of the surface density 
S for ass =0.1 (left) and ogg = 0.01 (right). The numbers are 
given in Table 1. Pm=l, B™''* = 1 
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Figure 5. The magnetic-disc height in comparison with 
nonmagnetic-disc height as a function of a) the vertical magnetic 



field B" 



/3 = 2 - 10"'^ (solid) and M 



4 - 10^ 



g/s (dot- 



ted), oss=0.01; b) the inclination angle i for ass=0-l (solid), 
oss=0.01 (dotted). Pm=l, B^°'^' = 1 



theless, this fact is only a rule with exceptions. Note that 
magnetic discs can be thinner than nonmagnetic discs for 
both large 5™'* and Pm. However, we should care the radial 
flow Mach number. The radial flow can become supersonic 
for H < Ho for large enough ass. This problem completely 
disappears for small ass. The disc thickness decreases with 
increasing radial magnetic field and vertical magnetic field 
due to magnetic stresses (see Fig. ^ . It is interesting to note 
that the accretion rate is getting larger for thinner disc. The 
disc thickness weakly depends on Pm as well as ass- 
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Figure 6. The inclination angle i of magnetic field lines to the 
rotation axis for _B^'"''=1. LEFT: toroidal field fixed (/3 = -9 ■ 
10""'), RIGHT: column density fixed S =1.4 g/cm^ 

4.3 Inclination angle i 

We find that for rather low /? the inclination angle i can 
be larger than the critical value of 30° (Blandford & Payne 
1982) for jet launching even for Pm of order unity. This 
result is in accordance with our calculation for polytropic 
magnetic disc (paper I). Moreover, fog given surface density 
the inclination angle hardly depends on Pm (Fig||) also in 
accordance with the results of paper I. This means that the 
radial velocity is proportionate to 1/Pm in this case. Never- 
theless, the inclination angle depends strong on Pm for given 
accretion rate. The inclination angle increases for decreasing 
ass for fixed accretion rate and decreases for fixed surface 
density. Such a behaviour is the result of the behaviours of 
accretion rate and surface density with changing ass- 

According to (^) and (pj]), the accretion rate is con- 
nected with /3. So, for given 5^°"^' there is a direct relation 
between the /3 and the resulting inclination angle i. Table 1 
present the numbers. Our model yields high accretion rates 
and high inclination angels already for rather small /3. The 
larger the fi is the higher the i. According to (|^), the /3 in- 
creased inwards. We can conclude that the Blandford-Payne 
condition is fulfilled most easily for the inner part of the ac- 
cretion disc. 



5 CONCLUSION 

The vertical structure of accretion discs with an imposed 
vertical magnetic field has been considered. The magnetic 
field energy is supposed as to be in balance with the accre- 
tion flow energy. The angular momentum transport is fully 
provided by the magnetic held, the Reynolds stress does here 
not play any role. 

The turbulent viscosity might only be important for the 
energy balance of the disc but the calculations demonstrated 
also the dominance of the Joule heating. The magnetic field 
is, therefore, the essential feature of the model and the equa- 
tions do not change to standard a-disc equations in the small 



magnetic field limit. Stehle & Spruit (2001) results is also 
confirmed the existence of magnetic field induced accretion. 
The results actually demonstrate the close relation between 
magnetic field dragging and the vertical structure of thin 
accretion disc. The interaction of the magnetic field with 
the disc can drastically change the disc structure as well as 
the configuration of the magnetic field. The angular velocity 
differs from the Keplerian one. This difference is, however, 
relatively small because the magnetic energy is small com- 
pared with the gravitational energy in our model. The radial 
velocity, however, can increase drastically becoming compa- 
rable to the sound speed for some models. This can lead to a 
strong amplification of the accretion rate for a given column 
density (Fig.^. 

The accretion rate is connected directly with the sur- 
face toroidal magnetic field in our model. Nevertheless, even 
small toroidal magnetic field at the surface (/? ^ 1) is 
enough to increase the accretion rate drastically. The accre- 
tion rate can be as amplified drastically (large i and small 
Pm), so suppressed (small i and large Pm). 

The magnetic disc is basically thicker than the nonmag- 
netic one for typical parameters values. Nevertheless, we can 
not exclude that the magnetic disc can became thinner than 
a nonmagnetic one for some parameters (e.g. large B™"^*). 

However, the most surprising results of our calculations 
concern the inclination angle i of the magnetic field lines 
to the rotation axis. We found that i) already rather small 
toroidal field component B^°^ can produce inclinations ex- 
ceeding the critical value of 30° and ii) this effect is almost 
independent of the magnetic Prandtl number for given sur- 
face density Fig. (^. 

In previous studies (Lubow, Papaloizou & Pringle 1994; 
Reyes-Ruiz & Stepinski 1996) large radial inclinations could 
only be obtained for Pm > 100. After our results it also 
holds for Pm=l. The difference is due to the fact that the 
previous studies neglect the magnetic field influence on disc 
structure and used radial velocity from standard accretion 
disc theory. 

For a given accretion rate, Af , the larger /3 lead to the 
larger i and to the smaller E (Fig. As after (^) the larger 
(negative) /3 exist in the innermost accretion disc region, 
we have there the smaller column density and the stronger 
radial inclination of the field lines. The jet launching should 
thus be concentrated to the inner region of an accretion disc. 

Only Kramers opacity was used in this paper but it 
might be very interesting to apply other opacity laws to re- 
formulate the stability problem of the structure of accretion 
discs. 
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